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Abstract 

Previously, we proposed a model of low energy Affleck-Dine leptogenesis in the 
context of thermal inflation. The lepton asymmetry is generated at the end of 
thermal inflation, which occurs at a relatively low energy scale with the Hubble 
parameter somewhere in the range 1 keV < H < 1 MeV. Thus Hubble damping 
will be ineffective in bringing the Affleck-Dine field into the lepton conserving 
region near the origin, leaving the possibility that the lepton number could be 
washed out. Previously, we suggested that preheating could damp the amplitude 
of the Affleck-Dine field allowing conservation of the lepton number. In this paper, 
we demonstrate numerically that preheating does efficiently damp the amplitude 
of the Affleck-Dine field and that the lepton number is conserved as the result. In 
addition to demonstrating a crucial aspect of our model, it also opens the more 
general possibility of low energy Affleck-Dine baryogenesis. 
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1 Introduction 



The origin of the observed baryon asymmetry is a fundamental question in modern 
cosmology. Forty years ago Sakharov showed that three conditions must be satisfied to 
generate the asymmetry [TJ, and since that time various models have been proposed to 
explain it. But we still don't have any idea which, if any, is correct, because baryogenesis 
usually depends on unknown physics beyond the reach of current experiments and on 
the unknown early history of the universe beyond the reach of current observations. 

A successful cosmological history must produce necessary relics such as the baryon 
asymmetry, but must also avoid or dilute unwanted cosmological relics such as graviti- 
nos [2] and moduli [3], which can destroy the successful predictions of Big Bang nucle- 
osynthesis jl] or over-dominate the universe depending on the scale of supersymmetry 
breaking. Thermal inflation (5j [6] provides a compelling solution to the unwanted relic 
problem. It occurs after the usual high energy inflation, at an energy scale between the 
intermediate and electroweak scales, releasing enough entropy at a low enough energy 
scale to dilute the unwanted relics. 

However, the usual baryogenesis scenarios do not work well in the context of thermal 
inflation because of the large entropy production at a low energy scale. Thus a new model 
of baryogenesis that fits naturally with thermal inflation is desirable [7J |8]. In a previous 
paper [9], we proposed a simple model based on the superpotential [j] 

W = \ u QH u u + X d QH d d + X e LH d e + \^ 2 H u H d + l -\ v (LH u f + ^A 4 (1) 

where is the flaton whose roll out from the origin ends thermal inflation and gener- 
ates the Minimal Supersymmetric Standard Model (MSSM) /x-term. The flaton must 
also have some unsuppressed couplings to the thermal bath, in order to hold it at the 
origin during the thermal inflation, but the extra fields involved in these couplings ac- 
quire intermediate scale masses after the thermal inflation and hence are only weakly 
constrained, so we do not include them explicitly. The baryogenesis is an Affleck-Dine 
(AD) [10] type of leptogenesis using L, H u and H d , and triggered by the roll out of 
the flaton. In Ref. [5], we demonstrated the generation of a lepton asymmetry using a 
simple homogeneous numerical simulation. In the usual high energy AD scenarios, the 
amplitude of the AD field would then be damped by the Hubble expansion, bringing 
the field into the symmetric region of the potential near the origin, ensuring that the 
generated lepton number was conserved. However, in our low energy AD scenario, the 
Hubble damping is negligible. Instead, we suggested preheating and thermal friction 
might damp the amplitude of the homogeneous mode, but testing this idea was beyond 
the scope of our homogeneous numerical simulation. 

During preheating, energy is transferred from a homogeneous oscillating mode to in- 
homogeneous modes, resulting in reduction of the amplitude of the homogeneous mode. 
Well known mechanisms of preheating include parametric resonance [UJ and tachyonic 
instability [12]. Usually, it is regarded as an intermediate process between the end 
of inflation and full reheating, i.e. thermalization. However, it can be relevant to the 

: We set c = h = 8ttG = 1 throughout this paper. 
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evolution of any condensate, in particular AD condensates. Especially, it may be a 
crucial effect for low energy AD scenarios, providing an alternative to the usual Hubble 
damping. 

In this paper, we perform three dimensional lattice simulations to numerically in- 
vestigate the effect of preheating on the amplitude of the AD fields, and hence on the 
conservation of the lepton number, in our model. In Section [2j we briefly review our 
model. In Section |3j we describe our scenario for the generation and conservation of the 
lepton number asymmetry. In Sections [4] and [5j we describe and present the results of 
our numerical simulations. In Section [6j we discuss interesting points, limitations of our 
current work, and work to be done in future. 

While this work was in progress, and after some of our initial results were reported 
[13) . another work on the same topic appeared [Hj. Their simulation was one dimen- 
sional, which restricts the number of modes available for preheating, and they chose 
different initial conditions and parameter^} This may explain why their results look 
somewhat different from ours, though they come to the same conclusion that the model 
seems to work. Unfortunately, there are insufficient details of their simulation in their 
paper for us to reproduce their results. We hope to clarify these issues after this paper 
is published. 



2 The Model 

Our model is based on the superpotential of Eq. ([I]). The fourth term generates the 
MSSM //-term when the flaton (f> settles down to its vacuum expectation value of the 
order of 10 10 to 10 12 GeV. The fifth term generates light Majorana masses for the 
neutrinos, and the sixth term stabilizes (p. We assume 

A M ~ A < A„ (2) 

to give the correct values for /i and the neutrino masses, assuming a flaton mass of the 
order of the electroweak scale. 

Initially, <ft is held at origin and so the //-term vanishes. Therefore H u and lack 
the usual MSSM //-term contribution to their mass-squareds. This allows the possibility 
that some supersymmetric flat directions involving H u or Ha, in particular LH U or 
H u Hd, may initially be unstable, although they are stable in our vacuum. Thus we 
consider that LH U , H u Ha and <fi may have nonzero values, and assume that the other 
fields remain at the origin [^J Then, truncating to a single generation for simplicity, gauge 
fixing and imposing the D-term constraints, we can parameterize those flat directions 
as 

Mo) ■ *-(£) ■ *-(?) < 3 > 



2 Their choice of parameters violates the MSSM constraint (rn%^ + |/i| 2 j (jn>% d + I MM < I^m| 2 ; but 
this is not important at this stage of simulations of our model. 

3 The non-trivial consistency of this assumption was checked in Ref. [9]. 
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with the remaining D-term constraint 



D=\h u \ 2 -\h d \ 2 -\l\ 2 = (4) 
and associated gauge constraint 

j = J2 rQD^-jDMQ^P = (5) 

ip£{h u ,h d ,l} 

which are not easily integrated out. The charges are Qh u = h u , Qh d = —h d and Ql = —I, 
and D a = d a — iQA a is the gauge covariant derivative. The gauge field A a is heavy and 
so not dynamical. 
The potential is 



V = V + m' Hu \KY + m^hdY + mi\l\* + m, 
+ 



A^\^4> 2 h u h d + ^A u \J 2 h 2 u + ^A^cj) 4 + c.c. 



+ \\f,(p 2 h d + XJ 2 h u \ +\X^(f) 2 h u \ +\XJh 2 u \ + |2A M (M„/i d + X^ (6) 

Vq is adjusted to give zero cosmological constant and all the other supersymmetry break- 
ing parameters are of the order of the electroweak scale. We denote a typical electroweak 
scale mass by m. 
We assume 

ml < (7) 
so that the thermal inflation ends by rolling away. We assume 

m 2 L + m 2 Hu <0 (8) 

so that LH U can roll away from the origin before the end of thermal inflation, providing 
the initial condition for our AD leptogenesis. This nontrivial assumption results in 
there being a deeper non-MSSM vacuum [JS] so that our MSSM vacuum is meta-stable. 
However, as explained in [9], this is phenomenologically and cosmologically consistent. 
For simplicity, we neglect any field dependent renormalization of the supersymmetry 
breaking parameters. 



2.1 Physical parameters of the potential 

In this section we give formulae for some physically important parameters of the poten- 
tial. Note that m\ < and m\ + m 2 Hu < 0. 

The magnitude of the flaton vacuum expectation value 0o is given by 



12m 2 + \A^\ 2 + {At, 
6|A^| 

The potential energy density at the origin is 



(10 10 to 10 12 GeV) 2 (9) 



H = ^(-"4+|A^| 2 )|0o| 2 ( 10 ) 



The flaton mass-squared eigenvalues at = O are 



m„ = 6 



2 \ - 2m 2 
o| zm </> 

2 1 2 , n o 



m a = 6|A^| + 2m 



(11) 
(12) 



where the subscripts r and a denote the radial and angular eigenvectors. 

When = hd = 0, the magnitudes of / and h u at the minima of their potential are 
given by 



— I^moI — 



-6 (m 2 L + m 2 H J + \A U \ 2 + \A V \ 
6|AJ 



(10 7 to 10 9 GeV) 2 (13) 



When settles down to (fi , the MSSM /i and B parameters are generated 



B 



and the AD mass-squared eigenvalues at the AD origin become 



m ln u 



1 I 1 2 | 1 / 2 i 2 \ 

1 



(14) 
(15) 

(16) 
(17) 



The mass-squared along the LH U flat direction changes its sign at the critical value C 
given by 



(Pc 



- (m 2 L + m 2 Hu ) 



To estimate the potential energy available in the AD sector, we use 



v 1 = v{^h u M )-v = 

V 2 EE V(<f) ,h u0 ,0,l ) 

where V((f>, h u , hd, I) is the full potential. 



(m 2 L + m 2 Hu ) + \XJl 



2|2 



\h\ +Vi 



(18) 

(19) 
(20) 



3 Leptogenesis 
3.1 Initial conditions 

Initially, all the fields are held at the origin due to the finite temperature effective 
potential. We assume lh u leaves the origin first. It will do so via a first order phase 
transition in which rare bubbles expand to Hubble volume size before percolating, since 
the expansion rate H ~ m \ <f> \ is much less than the scale of the phase transition m. As 
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lh u rolls away, it gives an increasing mass to the fields it couples to. These then decay 
and hence are eliminated from the thermal bath, but deposit their energy, including that 
of their extra mass gained from lh u , back into the thermal bath. This could lead to two 
possibilities for 0's dynamics, depending on 0's couplings and the detailed dynamics. 

Some fields are removed from the thermal bath, and if these fields couple to then 
0's finite temperature effective potential could be reduced, possibly triggering 0's phase 
transition. In this first case, lh u and roll away at the same time, though with lh u in 
the lead if |m|| < |m| + rn 2 Hu \ /2. 

However, the temperature of the thermal bath is raised which would have the op- 
posite effect, tending to prolong the thermal inflation. In this second case, lh u would 
have time to settle to its minimum, preheat and maybe even decay or be diluted by 
the remaining thermal inflation, before <fi undergoes its own first order phase transition, 
again with rare bubbles expanding to Hubble volume size before percolating. 

In this paper, we consider the second it seems simpler. The initial condition 

for the leptogenesis will then be 0's phase transition ending the thermal inflation, with 
lh u in some initial state near its minimum. In this paper, we will only attempt to 



crudely mimick 0's first order phase transition, as described in Section 4.2 A bound on 
the initial state of lh u can be estimated as follows. 

After lh u rolls away, the energy density in the AD sector is initially 

Padq ^ m 2 \l \ 2 (21) 

The AD sector will decay to the thermal bath at a rate 

Tad > f- 2 (22) 
I to | 

and the rate of transfer of energy from the AD sector to the thermal bath will be 

1 c?Pbath . Pad Tad , oq s 
+ 4p b ath = — 77 (23) 



Hq dt Hq 
Therefore the energy transferred in the first Hubble time 

» PadoFad v. m 4 

Ap bath £ > — (24) 

will be much larger than m 4 , and so the decay of the AD sector will be able to maintain 
the thermal inflation. The thermal inflation will end when pbath ~ fn A , i-e. after 

Pad < ^ < Pado ~ (10 2 ~ 10 8 ) m 4 (25) 



AD 



Thus, we can expect the initial fluctuations in lh u to be fairly small and well pre- 
heated if not thermalized. However, the leptogenesis seems not very sensitive to this 
initial state, see Section [5j 
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3.2 Generation of the lepton asymmetry 



When the temperature of the thermal bath drops to around the electroweak scale, it is 
(/>'s turn to roll away from the origin, ending thermal inflation. As rolls out, hd acquires 
linear terms in its potential which force it to become non-zero. When becomes of the 
order of its vacuum expectation value, h u acquires an additional contribution to its 
mass-squared from the term |A M 2 /i u | 2 , bringing it, and I and hd, back in towards the 
origin. At the same time, the phase of lh u is rotated by the term {\^4> 2 hd)* X u l 2 h u + c.c, 
generating the lepton asymmetry 

i*i - w 

n L = , (26) 

Note that the nonzero hd is necessary to generate the lepton asymmetry, and is also 
essential to stabilize the dangerous field directions Qd and Le which would lead one to 
the deeper non-MSSM vacuum 0. 



3.3 Conservation of the lepton asymmetry 

In the usual high energy AD scenario, the generated asymmetry is conserved because 
the Hubble expansion damps the kinetic energy of the AD field, bringing it into the 
symmetric region of the potential near the origin. However, in low energy AD scenarios 
the Hubble damping is negligible and the asymmetry could be washed out by rescattering 
off the asymmetric part of the potential away from the origin. Fortunately, the amplitude 
of the AD fields can also be damped by the transfer of energy to gradient energy and 
other fields. In this paper we focus on the former effect, which generally goes by the 
name of preheating [TTJ [T2"] . 

In our model, the gradient energy of the flaton initially grows very rapidly as the 
flaton rolls away from the origin, due to the tachyonic potential in the radial direction 
and the divergence of trajectories in the angular direction. This rapid growth flattens off 
as the flaton settles down to oscillating about its minimum, with modest and declining 
growth of gradient energy due to the non-quadratic terms in the potential around the 
minimum. The formation and decay of flaton top logical defects also has an effect. 

The AD potential has a strong dependence on the flaton, so inhomogeneities in the 
flaton field are transferred to the AD fields, leading to rapid preheating of the AD fields. 
In addition to this, preheating due to the couplings among the AD fields themselves is 
likely to be significant, but it is difficult to analyze analytically due to the complicated 
form of the AD potential. 

The energy of the AD fields can also be transferred to other fields. When an AD 
field passes near the origin, some fields which couple to it, such as quarks, may become 
sufficiently light that they can be thermally populated. As the AD field moves away from 
the origin, the fields become heavy again, extracting energy from the AD field. This 
energy is then transferred to the thermal bath, initially composed mostly of gluons, as 
the fields decay. This results in a friction like effect, which we call 'thermal friction', on 
the motion of the AD fields when they pass near the origin. The efficiency of thermal 
friction will depend on how close the AD field has to pass to the origin for the quark 
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or other field to become light enough to be thermally populated, and also on how high 
the temperature is. How close the AD field has to pass to the origin will depend on 
the couplings and the temperature of the thermal bath. The initial temperature of the 
thermal bath may be relatively small, perhaps electroweak scale, but the MSSM has 
a large hierarchy of Yukawa couplings. This should allow first generation quarks and 
leptons to be thermally populated when an AD field passes only moderately close to the 
origin. As the temperature rises, quarks and leptons of the higher generations will also 
be able to be populated. 

Acting against these effects, the flaton dependence of the AD potential will tend 
to transfer energy from the flaton to the AD fields. The initial flaton energy density 
Ptp ~ m 2 1 0o 1 2 is much larger than the initial AD energy density pad ~ m2 | | 2 , so there 
is potentially a lot of energy to be transferred. However, one would expect the energy 
flow to decrease as the amplitude of the flaton decreases due to preheating, energy 
transfer and eventually Hubble damping. 

The efficiency of these effects will be crucial in determining whether our lepton 
number asymmetry survives. However, it is unclear how efficient any of these effects is 
and even which of them is dominant. So, to start somewhere, in this paper we analyze 
the preheating effect numerically, which also includes the energy flow from the flaton to 
the AD fields. 



3.4 Estimation of the final baryon number asymmetry 

Once the amplitude of the AD fields has been damped sufficiently for the lepton vio- 
lating terms \A u \J 2 h 2 u and (X^hd)* \J 2 h u + c.c. to become negligible, B — L will be 
conserved. We denote this initial time by t\. The AD fields will complete their decay 
to the thermal bath while the temperature is still above the electroweak scale, and so 
sphaleron processes will convert our lepton asymmetry to a baryon asymmetry. Finally, 
the flaton decay completes, diluting the baryon asymmetry. We denote this final time 
by t{. 

The final baryon asymmetry is given by 

n B n B _ L (t { ) n B - L {ti)T(t { ) n L (ti)n B - L (t { )T(t f ) 



s s(t { ) p(t f ) n B _ L (£i) p(tf) 

'm LHu n L (ti)\ fV 2 \ ( T(t f ) \ ( n B - L (t { ) V 



(27) 



V 2 J \V J \m LHu J \n B - L (U) p(t { ) 
where p and T are energy density and temperature of the universe. We use the factor 

m LHu n L (t^ 



Vo 



(28) 



to give a dimensionless measure of the efficiency of the lepton asymmetry generation 
and preservation. It can be interpreted roughly as the ratio of lepton number density 
to the total AD number density when the AD fields start oscillating. We expect 

^ ~ 10- 10 to 1(T 2 (29) 
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and 

~ 1(T 3 to KT 1 (30) 

™lh u 



as described in Section 2.1 and Ref. [9]. The last factor in Eq. (27) depends on the 
equation of state of the flaton. If the preheating of the flaton only populates modes 
with k ~ m, then the flaton would have p/p <C 1. The energy flow from the flaton 
sector to the AD sector could also affect this factor, but, although this flow is likely 
significant from the AD point of view, we guess it is not so significant from the flaton 
point of view. With these assumptions, one expects 

ns-L{t{) V 



n B -L(ti) p(t f ) 



1 (31) 



Thus, we obtain 



^ ~ (10"" to 10- 3 ) (32 ) 
s V 2 



compared with the observational result 



^ ~ KT 10 (33) 
s 



4 Numerical Simulation 

We performed a three dimensional lattice simulation with periodic boundary conditions 
using the algorithm described in the Appendix. The potential was that of Eq. ^ with 
constraints Eqs. Q and ([5]), and canonical gauge invariant kinetic terms. The constraint 
Eq. Q was modified to 

D = e 2 (34) 

to cut off the singularity at h u = 0. The physical parameters and fields were rescaled by 
a typical electroweak scale mass m and a typical flaton or AD field expectation value 
Mti or Mad, as described in Ref. [9]. 



4.1 Simulation parameters 

The lattice volume L 3 , number of lattice points iV 3 , time step At, and D-term constraint 
singularity cutoff e, were taken as 

L = 50m~ 1 , iV = 128 , At = 4 x Kr 3 ™" 1 , e = 5 x 10~ 3 M AD (35) 

We tested our results using different values of these numerical parameters. Limited 
computing power constrained us to iV < 128. 

The fc-modes allowed by this lattice are k = \/k\ + k\ + k\ with 
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where n$ = 0, 1, . . . , N/2 with rij = 0, N/2 having degeneracy 1 and the rest degeneracy 
2. This spans the range 



0.13m = — < k < ^— — = 14m 

Lj Lj 

The physical parameters were taken as follows. 



(37) 



ml = -0.217m 2 



m H 



-1.510 m 2 



m 2 H = 3.533 m 2 



m 2 r = 1.223 m 2 



(38) 

The mass-squareds were given with values up to three decimal places to avoid any 
accidental resonances and also to fit the MSSM constraints described in Ref. [§]. 



\A„ 



1.0 m 



\A, 



I A, 



l.OmM- 2 



IA, 



0.5 m 



2.0 mM^ 



\AJ = 0.25 m 



|A„| = 0.2 mMll 



M AD = 10" 2 M TI 
arg {A^) = arg (A^) = arg (A v \ v ) = 



(39) 

(40) 
(41) 
(42) 



Note that the phases of Eq. (42) are not physical and can be adjusted by field rotations. 
The CP phases were taken as 

aig{-A*A„) = 



71 
71 
71 



71 

20 

71 

20 



CP+ 
CP0 
CP- 



arg(-A;A M ) = 



For the choice of parameters above, the formulae of Section 2.1 become 



0.70 M T i 



m 2 = 6.0 m 2 



1.7 x 10" 2 |0 O | 
= 0.97m 



Vq = 0.23 m 2 



m 2 = 5.2 m 2 



\B\ 



= O.74|0 O | 
1.5 m 



m 



LH U 



0.32 m 2 



ni 



2 rt tt = 3.4 m 2 , m'rr 



0Mm z 



Vx = -2.6 x 10" 4 V 



V 2 = 8.7 x 10" 4 V Q 



(43) 
(44) 

(45) 

(46) 
(47) 
(48) 
(49) 
(50) 
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4.2 Initial conditions 



The real initial conditions for the leptogenesis includes the realization of the first order 
phase transitions of lh u and 0, or at least the first order phase transition of and the 



state of lh u at that time, as discussed in Section |3.1[ During the transitions, bubbles 
of non-zero field value nucleate, expand and percolate. The typical size of a bubble 
just before percolation is expected to be slightly smaller than the Hubble radius. This 
introduces a scale much larger than the inverse electroweak scale typical of the rest of the 
dynamics, making it difficult to simulate both the phase transitions and the leptogenesis 
dynamics simultaneously. We leave detailed analytic and numerical investigation of these 
initial conditions to a subsequent paper, and here just use a simple ansatz that covers 
a range of possibilities, hopefully without missing too much relevant physics. 

We do not include a thermal bath as after the fields roll away from the origin, they 
decouple from the thermal bath, apart from effects like thermal friction which we plan 
to consider in a subsequent paper. 

We set the initial conditions as 

= A0 + 50 (51) 

1 = k + 51 (52) 
h d = 5h d (53) 

with the initial conditions for h u determined by 

D = e 2 (54) 

D = (55) 

arg/i u = 7r/4 (56) 

A) = (57) 

We choose the initial homogeneous displacement of to be either A0 = or A0 = 4m. 
A0 = corresponds roughly to a second order phase transition, and provides one 
extreme. A finite value of A0 is meant to correspond roughly to the displaced field 
inside a bubble. D is given by Eq. Q and A = leads via Eq. (|5p to 

= kk - KK _ h* d h d - h* d h d _ in - in = Q 

i i i 

5(p, 51 and 5h d are inhomogeneous fluctuations taken to be the combination of vacuum 



and thermal fluctuations at a temperature T ~ m. Since, as discussed in Section 3.1 
in some cases SI may be quite large, we also test the effects of an initial spectrum of 51 
fluctuations with a thermal spectrum of temperature T ~ 10 6 m up to the lattice cutoff 



fc mH ~ 10m, roughly corresponding to a preheated state with \l (\5l\ 2 ) ~ 10 2 \l 
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5 Results 



5.1 Lepton number asymmetry 

Our most important results are Figs. [2] and [I] showing the lepton number asymmetry 
plotted as a function of time. In Figs. [2] and [4j the red and blue lines show that a 
lepton number asymmetry is generated, with sign depending on the sign of the CP 
violating phases. The green line is the result of a simulation with CP conserving phases 
used to estimate the noise level due to the random initial conditions and finite size of 
the lattice. In the case of A</> = 4m, the angle of initial displacement can contribute 
to the CP violation and hence affect the lepton asymmetry, as shown in Fig. [3j We 
average over this initial angle, corresponding to summing over bubbles with different 
initial flaton phases, to obtain the lepton asymmetry of Fig. [2] A common feature of 
Figs. [2] and [4] is that, although some of the initially generated asymmetry is washed out, 
a significant fraction is conserved to late times. 

The initial conditions, A0 = and A0 = 4m, correspond to a range of initial A<p 
so these results suggest that conservation can be achieved for a wide range of initial 
conditions. Moreover, Fig. [5] shows that the lepton asymmetry is not sensitive to the 
state of lh u at the time when the flaton starts to roll away. This suggests that the 
fluctuations of the AD fields acquired from those of the flaton through their couplings 
dominate over any initial fluctuations of lh u . 

The lepton asymmetry declines slowly at late times. This is due to the preheating 
becoming saturated and energy flow from the flaton sector, as will be discussed in 
Sections 15.21 and 15.31 

In Figs. [4] and |5j the CP+ lepton asymmetries are enhanced after a short period 
of settling down. It is unclear what this is due to, but it does not seem to be harmful 
because it disappears leaving a large enough lepton number asymmetry. The spatial 
distribution of the lepton number density shown in Fig. [7] shows a large ball type struc- 
ture, which seems to be responsible for the strange behavior in Fig. |4} We suspect it 
may be a quasi-stable Q-ball [HI [13 [T5] , but it is unclear. We postpone proper analysis 
of this to future work. 

Thus we have demonstrated the central aim of this paper, that preheating can damp 
the amplitude of the Affleck-Dine fields sufficiently to conserve the generated asymmetry. 
We hope to extend this result to more realistic initial conditions in a future paper. 

5.2 Preheating 

The dispersions and spectra of the fields give a more detailed picture of the dynamics 
leading to the conservation of the lepton asymmetry. The field amplitudes in Figs. [8] and 
[9] (top left of each figure) show the AD fields settling down into the symmetric region 
of the potential near the origin, which leads to the conservation of the lepton number. 
This is a consequence of very efficient preheating as shown in the graphs of the field 
dispersions (left middle and bottom) and mean squares about vacuum expectation value 
(top right) in the figures. While the decreasing behavior around t = 100 in the graphs of 
mean square about vacuum expectation value is due to energy flow from homogeneous 
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modes to inhomogeneous modes, the similar behavior in the dispersion graphs is due to 
energy flow to higher /c-modes. 

As is clear from Figs. [2] and |4| and the graphs in Figs. [8] and |9j the two cases of 
A(p = 4m and A0 = have different initial transient behaviors. This is due to the 
difference in the efficiency of the preheating of the flaton. For A0 = 4m, the flaton has 
an initial homogeneous displacement larger than the fluctuations, so that it rolls away 
along a specific direction with relatively small fluctuations, leading to slower growth 
of gradient energy. This results in the flaton oscillating several times before settling 
down to its minima, and causes the violent initial behavior of the AD fields due to 
their flaton-dependent mass-squared's whose sign also depends on the flaton field value. 
However, for A<ft = 0, the flaton rolls away from the origin with initial fluctuations in all 
directions. This leads to rapid growth of gradient energy, resulting in the flaton settling 
down to its minima within a single oscillation. 

Note that the large initial flaton dispersion and its drop around t = 100 in Fig. [9] is 
due to the formation and annihilation of flaton strings and walls, which can be seen in 



Fig. 11 In contrast, Fig. 10 shows that the A0 = 4m case is initially dominated by a 
single domain, as would be expected due to its homogeneous displacement. Our model 
has a Z 4 gauge symmetry identifying the four flaton vacua [9J. However, a discrete gauge 
symmetry is more difficult to simulate than a global symmetry, so in our simulation we 
used a Z4 global symmetry, and relied on the random initial conditions and finite size 
of the lattice to get the strings and walls to annihilate, in accord with the dynamics 
expected for the Z 4 gauge symmetry. 

After the initial transient stage, the amplitudes of the AD fields start to grow again. 
This is because the energy flow from the flaton sector becomes greater than the rate 
of energy dispersion to higher /c-modes due to preheating and thermalization. This is 
the reason why the lepton asymmetry in Figs. [2] and [4] slowly decreases at late times. 
The difference in the late time rate of increase of the AD amplitudes in Figs. [8] and [9] 
is due to difference in the rate of energy flow into the AD sector, which is caused by 
the difference in field amplitudes at the time when preheating is effectively shut down. 
So the efficiency of preheating at the initial transient stage determines the late time 
behavior of the graphs. This is shown well by the lattice cutoff dependence shown in 



Fig. 17 



The spectra in Fig. 12 show the amplification of fc-modes through preheating, which 
damps the homogeneous motion of the flaton and AD fields. In the flaton sector, we 
initially have tachyonic growth of low fc-modes as the flaton rolls away from the origin. 
After this initial stage of preheating due to instability, the flaton oscillates around its 
minimum. At this stage, nonlinear effects populate slightly higher £;-modes of the flaton. 
At late times, the flaton spectrum is cut off at high k, which is consistent with theoretical 
expectations [12]. The preheating of the flaton sector is transferred to the AD sector 
through the flaton-dependent potential of the AD fields. Remarkably, AD fc-modes 
higher than the physical cutoff in the flaton sector are also rapidly populated hitting the 
high k lattice cutoff before t = 100 with rapid energy flow from lower to higher fc-modes. 
As a result, the AD lattice modes are fully thermalized for most of the simulation. 

Once all the lattice modes have been thermalized, preheating is truncated artificially 
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in the simulation, even though it may continue in reality. This suggests that the relia- 
bility of the lattice results will be compromised at late times, as energy that physically 
should be flowing into the higher k modes is artificially, numerically shunted into the 
modes that are present in the grid. This effect should artifically weaken preheating and 
we believe it accounts for some of the late time decrease in the lepton asymmetry seen 
in Figs. [2] and [4} 

We have performed experiments with coarser lattices (N = 64, 32) and shown that 
preheating in the AD sector becomes stronger as the UV resolution is increased, con- 
firming our expectation that our limited resolution is artificially damping the effects 
of preheating. (These experiments also show that the grid resolution is not affecting 
the flaton sector.) Our results thus present a lower bound on the effects of preheating. 
Thus in reality one should expect more preheating in the AD sector and hence better 
conserved lepton number than the simulation shows. We postpone further study of this 
issue to future work. 



5.3 Energy flow 



Figs. 13 to 14 show the components of the energy of the flaton and AD fields as a 
function of time with total energy densities defined as 



PAD = { \D Q ll>\ 2 +Y,\Dd>\ 2 + VAD 



(59) 
(60) 



i=i 



where (...) denotes the lattice average and 



V FL (0) 

Vad(0,VO 



Vo + mi\ 

v-v FL 



+ c.c. 



3 I 2 



(61) 

(62) 



As <fi rolls away from the origin, its potential energy Vq is rapidly converted into 
kinetic and gradient energy. Even at late times, the flaton's potential energy is slowly 
being converted into kinetic and gradient energy, due to continued preheating. The 
flaton dependence of the AD potential initially generates potential energy for the AD 
fields, which is then rapidly converted to kinetic and gradient energy. At late times, 
the AD lattice modes have thermalized and the energy from the flaton sector flows into 
kinetic, gradient and potential energies of the AD fields. The total energy flow from 



flaton to AD fields is shown in Figs. 15 and 16 



The qualitative features of Figs. [15] and [16] can be understood from the energy flow 
rate into the AD sector 

(63) 

which is obtained by use of the equation of motion of <fi. We can separate the energy flow 
into three stages. At the first stage, the AD potential is lifted up causing the AD fields 
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to oscillate with large amplitude. Then, preheating spreads energy into higher /c-modes, 
reducing the field amplitudes. This results in decreasing energy flow rate. Finally, the 
AD sector preheating is saturated due to the lattice cutoff, so energy from the flaton 
sector starts being accumulated in a fixed number of modes. This causes increasing field 
amplitude and therefore increasing energy flow rate. Thus, the late time behavior of the 
energy flow depends on the efficiency of preheating, which in turn depends on the initial 
conditions and the lattice cutoff. The dependence on the initial conditions explains the 
different rates of decline in Figs. [2] and |4} 

We expect a smaller energy flow for a finer lattice, as shown in Fig. 17, but this will 
be the case only until a physical cutoff of the AD preheating appears. If the preheating 
becomes saturated due to a physical cutoff and the energy flow is still significant, then 
we may need the help of thermal friction which could provide another channel for energy 
to flow away. We plan to study these issues in a future paper. 



6 Conclusion 

Recently, we proposed a minimal extension of the MSSM which naturally gives rise to 
thermal inflation followed by AD leptogenesis. The expansion rate H at the end of 
thermal inflation is too low for the usual AD Hubble damping mechanism to conserve 
the generated lepton asymmetry, but we proposed that preheating and thermal friction 
may perform the same function. However, those processes are very complicated and not 
easy to study analytically even in a simple case. It is even difficult to know which of 
these processes will dominate since it depends on the details of the dynamics. 

In this paper we numerically analyzed the effect of preheating on the conservation 
of the lepton number asymmetry in our model. The lepton number asymmetry is gen- 
erated as the flaton rolls away at the end of thermal inflation. It will be conserved if 
the preheating of the AD sector is efficient enough that the homogeneous oscillation 
amplitudes of the AD fields are quickly reduced, bringing the fields into the symmetric 
region of the potential near the origin. 

The preheating process starts in the flaton sector by the magnification of the initial 
fluctuations, that include vacuum and thermal contributions, as the flaton rolls away 
from the origin and oscillates about the minimum of its potential. Since the flaton 
potential has negative curvature near the origin and somewhat larger angle-dependent 
curvature along angular directions far from the origin, very rapid tachyonic preheating is 
dominant initially. When the flaton oscillates around its minimum preheating continues 
due to the anharmonic terms in its potential (weak parametric resonance). 

The AD fields in our model have flaton-dependent masses, so their mass squareds 
vary with time and even change sign. This violently transfers inhomogeneities in the 
flaton sector to the AD sector. The preheating amplifies high fc-modes in the AD sector 
much more rapidly than in the flaton sector, so that the AD fields quickly thermalize 
over the range of momenta present in our lattice. If we were able to simulate with higher 
resolution we would expect to see even more efficient preheating. 

The flaton sector also feeds energy into the AD sector. The energy flow rate depends 
on the amplitudes of the fields and the time derivative of the flaton. This means that 
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if the flaton and AD fields are efficiently preheated, resulting in small field values, the 
energy flow rate will be small. Since the AD sector preheating is very rapid initially, 
the energy flow in the initial transient stage is not harmful in terms of lepton number 
conservation. But the preheating is quickly saturated, causing the energy flowing from 
the flaton sector to accumulate in a fixed number of modes. This causes an increase in 
the oscillation amplitude of the AD fields, which can wash out the asymmetry, as well 
as increase the energy flow rate. In this sense, the balance between the preheating in 
the AD sector and the energy flow from the flaton sector is crucial for the late time 
evolution of the lepton asymmetry. 

There is a large structure in the spatial distribution of the lepton asymmetry, which 
may be a quasi-stable Q-ball, but it is not harmful since it soon disappeared without 
much effect on the asymmetry. The origin of this object is not clear. Including the 
renormalisation of the AD masses should make such objects less stable, assuming it is 
a type of Q-ball. 

The outcome of all these effects is, as the simulations show, that the lepton asymme- 
try in our model is reasonably well conserved at late times due to the efficient preheating, 
and is not affected much by the initial state of the AD fields. The lepton asymmetry 
decreases slowly at late times but the rate is small enough for one to expect that suf- 
ficient lepton asymmetry will be preserved for long enough for other physical processes 
like thermal friction to start working and stop the decrease. 

There are several limitations to this work, which may affect the conservation of 
the lepton asymmetry. Firstly, all of our results may depend crucially on the initial 
conditions. We found that the initial state of the AD fields has little effect on the lepton 
asymmetry. The initial state of the flaton, however, affects the preheating of the AD 
fields and the energy flow, and thus may affect the conservation of the lepton asymmetry. 
A proper treatment of this issue will thus require simulating flaton bubbles. Once this 
has been done, it becomes possible to determine the regions of parameter space in which 
our model works. 

The limited size and resolution of the lattice also affect the preheating efficiency. 
Since the lattice cutoff is reached early in our current simulation, it restricts the pre- 
heating artificially. A larger and finer lattice would provide more available modes for 
preheating and thermalization such that it would weaken or even stop the decreasing 
behavior of the lepton asymmetry. 

Finally, we may have to consider thermal friction as another source of damping 
for the homogeneous mode, if it turns out that preheating is not efficient enough to 
conserve the lepton asymmetry by itself given realistic initial conditions for the flaton. 
We postpone all of these issues to future work. 

In summary, our simulation revealed that preheating is efficient and rapid enough 
to conserve the lepton asymmetry in our low energy AD leptogenesis model, replacing 
the Hubble damping of the usual high energy AD scenario. 
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A Appendix 

A.l Gauge invariant adaptive leapfrog with constraints 

In this appendix, we describe our algorithm for numerical integration of our dynamics. 
It is a gauge invariant adaptive leapfrog type algorithm with constraints. It is derived 
from a discrete action, and so is a variational integrator [Tj5] , which means that it exactly 
conserves the constraints and charges, and has good long term energy conservation. The 
adaptivity is particularly important to help deal with the singular configuration manifold 
induced by the D-term constraint of Eq. Q . 



A. 1.1 Lattice gauge invariance 

The remaining gauge field in our model, corresponding to the D-term and gauge con- 
straints of Eqs. Q and ([5]), is heavy and so not dynamical, but, unlike the other gauge 
fields, can not be easily integrated out. 

A gauge field is a connection relating variables at two different points in space-time, 
and so should be considered as a link variable living on a link joining two points. We 
will follow the usual way of defining covariant derivatives on a lattice in terms of link 
variables. 

Let's denote nodes by n and links by /. A link I connects from the node I— to the 
node /+. The links n— flow into, and the links n+ flow out of, the node n. The length 
of a link I is A/. A link tl is in the time direction and a link xl is in the space direction. 



Covariant derivative 

Ui = e lQAlAl (64) 

(2X&)j = 1 \ 1 (65) 

(Z>$)J(X>$), = — ^ + A2 l+ (66) 
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At 



Ax 



Figure 1: Our spacetime lattice. Hollow nodes are added adaptively to the time links. 
The kinetic terms are applied on the adaptive time links, the gradient terms on the 
space links, the potential on the solid nodes, and the D-term constraint on all nodes. 

where Q is the charge operator, A\ is the gauge field on the link I, and $ n is the complex 
scalar field vector at the node n. 



Gauge transformation 



-4 ; 



^1+ — \- 



A ' + A, 

e iQ{Xi+-M-)jj 

AQ(\i++Xi-)/2 



(67) 

(68) 

(69) 
(70) 



where A ra is the gauge transformation parameter at the node n. 



Gauge fixing We set 

by taking X a+ = X a _ - A U A U . 



Aa = 



(71) 



A. 1.2 Discrete action 

Our discrete action is the sum of kinetic terms on time links, gradient terms on space 
links, and potential and constraint terms on nodes, with spacetime volume weightings 
fl appropriate to their domains 

s = E n« - e (d$) x i ^ - E v <& - E (P n - e2 ) 

tl xl n n 

(72) 

where 

D n = $t Q$ n (73) 

We choose a regular spacetime lattice with spatial lattice spacing Aa; = L/N and 
time step At for the gradient and potential terms. This regular lattice is denoted by 
the solid nodes and the links joining them in Fig. [TJ Extra nodes, the hollow nodes in 
Fig. [TJ are added adaptively to the regular lattice's time links, and we use this adaptive 
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lattice for the kinetic and constraint terms. This helps to follow the dynamics on our 
constraint manifold which is singular in the limit e — > 0. 
For this lattice, the weightings are 



Q tl = (Ax) 3 A tl 
Vt xl = (AxfAt 
nl = 6. n (Ax) 3 At 



where 



A. 1.3 Constraints 



1 for n 
for n 



(74) 
(75) 
(76) 
(77) 

(78) 



Varying our action of Eq. (72) with respect to the Lagrange multiplier A n gives the 
D-term constraint 



Therefore 

= {dD\ 



= l -— = D - e 2 = -e 2 

A+ - A- ®] + Q$i+ - $\-Q$i- 



A, 



A, 



i \t 



(79) 



(80) 



(81) 



Varying with respect to the gauge field A\ gives the gauge constraint 
9u SS T $\_Ur x Q$i + - < 



n t 6 A 



Jl 



1A1 

I \t 



(82) 



2i 



Combining Eqs. (81) and (83) gives 



0=(dD + iJ)i = \U l ^+ + ^yQ(X?$) ; 



2$[ + f/ i 5 g(P$) i - (P$)[Q(P$) i A i 

2$J_e/~*q(x>$), + (z>$)Jq(p$)a 



(83) 

(84) 

(85) 
(86) 
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Charge spectrum with Q 3 = Q For $ = ($ + , $°, $ ) with 



g$+ = $+ , g$° = o , Q$- 



Eq. (82) simplifies to 



-t, 



Choosing the positive solution consistent with the limit Ai — > 0, we get 



Note that in our case of $ + = /i u and $ = (/i^, /), Eq. (79) gives 



We use Eq. (89) to determine U x i, and from Eq. (71) we have 

U tl = l 

A. 1.4 Dynamics 



Varying our action of Eq. (72) with respect to <&t gives 
1 5S 







where 



(Ax) 3 AU$t 



<W A 8D 

(V?$) - 5 (V 2 $) + 5 — + -5 A — 



At 



the hat indicates improper normalization due to the use of At instead of A tn ±, 

u:} + (v$) xn+ -ui_{v$) xn „ 



xn+ 



Ax 

- 2$ n + U xn - <$> xn . 



Solving for A n we get 



t 

(Ax) 3 At AnVn 



where 



(Ax) 2 



— — — c)V 
-n{(V^) n + (Wit (!#&)„ - 5. n nt — » 



™n = — r = Q$n 
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Therefore, defining 



and 



nn 



Eq. (92) becomes 



A. 1.5 Algorithm 



Pn = l~N n 

P n (pf&) n - 5. n P n {V 2 x §) n + 6 m P, 



dV n 



Eq. (86) gives 



N n U tn l(V<S>) tn+ = - 1 - [(P$)tQ(P$ 



tn+ V n \ 



n 



n<n/ n 



A 



and Eq. ( 100 ) gives 



tn+ 



P n U tn l(V$) tn+ = P n Ul_{V®) tn „ + 6 m P n (T>Z$) n At - 5. n P n -^ At 



Combining these gives 

u; n i(v$) tn+ = p n 



i QV 

Ul_(m>)i^ + S. n (V 2 x $) n At - 5. n — j At 



- - [(P$)tQ(P$)] 



tn+ \ n ] 



n 



rvn/ n 



A 



tn+ 



Implementing our gauge choice of Eqs. (71) and (|9T|) and defining 



7I~n — Pn 



dV n 



tn _ + 5. n (V 2 x $) n At - <W~4 At 



and 



Atn+ 



A r 



reduces Eq. (103) to 



tn+ — K n — — Hf n ^.n n A n 



while Eq. (65) reduces to 



$ tn++ = $„ + $tn+A 



tn+ 
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Eqs. (108) and (109) are our basic algorithm, but we still need to determine II tn+ . 



Contracting Eq. (108) around Q gives 



IW = (7rtQ7r) n - l -Ii tn+ [(^Qn) n + c.c] A n + hl 2 tn+ (n ] Qn) n K 2 n (110) 
Choosing the solution which is well behaved in the limit A n — ► 0, we get 



n 



tn+ 



and 



1 + 2 ( nkQiTn + C.C. ) A n + \ j 1 + 2 ( nlrQ-rin + c.c. ) A n - ■K J n Q'K n n 1 n ,Qn n Al 

(111) 

(112) 
(113) 



For Eq. (Ill) to remain well behaved we require 

1 



iv n Qix n + c.c. A n < 1 



\ 7r nQ 7r n n nQ n n\ A 2 < 1 

which are the conditions we use for the adaptive part of the algorithm. 



For a charge spectrum Q 3 = Q, Eq. (79) gives 



rJn > n^Qn = e 2 



Therefore a sufficient condition for Eqs. (112) and (113) is 



Kn| A tn+ < e 
but it is not necessary as typically rvn ^> e 2 . 

A. 1.6 Lepton number 

The lepton current is 



(114) 
(115) 



J 



& l+ u?L{v®) l -{v<s>)\u;>m + 



<$>\_U^L<$> 1+ - tf+UiLfy. 
z'Ai 



(116) 

(117) 
(118) 



where L/i u = 0, Lhd = and LI = I. For our gauge choice of Eqs. (71) and (91) the 
lepton density reduces to 



J, 



(119) 
(120) 
(121) 
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Our algorithm gives 



At 



J L - J L 

" xn+ " xn- 

Ax 



S.n + i 



(122) 



Therefore our algorithm exactly conserves lepton number apart from the lepton number 
violating terms in the potential. 



A. 1.7 Summary 



Our algorithm of Eqs. (108) and (109), with n n given by Eq. (105), (V^) n given by 



Eq. (95), U xn ± determined by Eq. (89) and TL tn+ given by Eq. ([Ill]), exactly conserves 
the D-term constraint 

D tl+ = D a - (123) 



the gauge constraint 



Ji = 



(124) 



and lepton number apart from the lepton number violating terms in the potential, and 
has good long term energy conservation, as is expected for a variational integrator [19J. 

It is also adaptive with respect to the kinetic and constraint terms, allowing efficient 
integration on our constraint manifold, which is singular in the limit that the cutoff 
e -> 0. 
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Figure 3: Lepton number density, averaged over the lattice, as a function of the initial 
phase of the flaton, for A(f> = Am. The Z4 symmetry of our model means that the range 
of initial angles is [0, 7r/2). 
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Figure 4: Lepton number density, averaged over the lattice, as a function of time, for 
A0 = 0. 
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Figure 5: Lepton number density, averaged over the lattice, as a function of time, 
for A<p = 0, CP+ and CP—, with different initial effective temperatures for lh u . We 
set an initial thermal spectrum up to the lattice cutoff k max ~ 10m, so that an effec- 
tive temperature of T ~ 10 6 m would correspond to a preheated state with fluctuation 
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Figure 6: Lepton number density as a function of the 2D lattice for time slices at 
t = 40, 60, 80, 100, 200, 500, for A0 = 4m, CP+, initial flaton phase 9$ = vr/4. 
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Figure 7: Lepton number density as a function of the 2D lattice for time slices at 
t = 40, 60, 80, 100, 200, 500, for A0 = 0, CP+. 
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Figure 8: Dispersions and variances averaged over the lattice, as a function of time, for 
A(j) = 4m, 9^ = 7r/4, CP+. Left top - field amplitudes in log scale; Left middle and 
bottom - dispersions in log and linear scales; Right top - mean square with respect to 
vev in log scale; Right middle and bottom - mean squares with respect to origin in log 
and linear scales. 
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Figure 9: Dispersions and variances averaged over the lattice, as a function of time, 
for A0 = 0, CP+. Left top - field amplitudes in log scale; Left middle and bottom - 
dispersions in log and linear scales; Right top - mean square with respect to vev in log 
scale; Right middle and bottom - mean squares with respect to origin in log and linear 
scales. 
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Figure 12: Kinetic energy spectra of flaton and AD fields, at i = 
0, 20, 40, 60, 80, 100, 200, 300, 400, 500 for CP+. Thermally normalized so that a thermal 
spectrum is flat and energy normalized so that the area under the graph is the total 
energy from top to bottom. Upper half: A0 = 4m, 9^ = 7r/4. Lower half: A0 = 0. 
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Figure 13: Energy densities of the flaton as a function of time, for A<p = 4m, CP+, 
00 = 7r/4. From top to bottom at t — 500: kinetic, gradient, potential. 
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Figure 14: Energy densities of the flaton and AD fields as a function of time, for A0 = 0, 
CP+. From top to bottom at t = 500: kinetic, gradient, potential. 
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Figure 15: Energy flow between flaton and AD fields, for A0 = 4m, CP+, 9$ = tt/4. 
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Figure 16: Energy flow between flaton and AD fields, for A0 = 0, CP+. 
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Figure 17: Dependence on lattice with L = 50 and iV = 32, 64, 128 for A0 = 0, 
CP+. Upper half: thermally normalized and energy normalized spectra in order at 
t = 0,20,40,60,80,100,200,300,400,500. Lower half: variances, total energies and 
lepton numbers densities , averaged over the lattice, as a function of time. The graph 
of the lepton number densities is including the cases of CP0 and CP— as well. 

38 



